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I. INTRODUCTION 

The idea of Gamow, proposed in the 1930s and 
1940s |jj about the production of elements through ther- 
monuclear reactions in the hot ylem during the early 
stages of the Universe expansion, has been developed dur- 
ing the last 60 years into an elegant famous theory of cos- 
mological nucleosynthesis, explaining quantitatively the 
inferred from observational data primordial abundances 
of the light elements Q . Thanks to that good accordance 
between theory predictions and the observational facts, 
we nowadays believe to have understood well the physical 
conditions of the nucleosynthesis epoch. Still, the uncer- 
tainties of the primordial abundances values extracted 
from observations yet leave a room for physics beyond 
the standard model. 

In this article we present a modification of the stan- 
dard model of cosmological nucleosynthesis (CN) - CN 
with neutrino oscillations. Our aim is twofold: (1) to 
construct a modification of CN using a more precise ki- 
netic approach to the problem of nonequilibrium neutrino 
oscillations and to illustrate the importance of such an 
exact approach, and (2) to determine the cosmologically 
allowed range for oscillation parameters from an accurate 
study of the oscillations effect on the primordial produc- 
tion of hclium-4, thus helping clarify the mixing patterns 
of neutrinos. 

The theme of neutrino oscillations is with us almost 
forty years, since the hypothesis for them was proposed 
by Pontecorvo || . They were studied experimentally and 
theoretically and their cosmological and astrophysical ef- 
fects have been considered in numerous publications Q 
as far as their study helps to go deeper into the secrets of 
neutrino physics and neutrino mass pattern. Nowadays 
there are three main experimental indications that neutri- 
nos oscillate, namely: the solar neutrino deficit Q (an in- 
direct indication), the atmospheric neutrino anomaly || 
(an indirect indication) and the LSND experiment re- 
sults (a direct indication). 

(a) Solar neutrino deficit: Already four experiments us- 
ing different techniques have detected electron neutrinos 
from the Sun, at a level significantly lower than the pre- 
dicted on the basis of the Standard Solar Model and the 
Standard Electroweak Theory. Moreover, there exists in- 
compatibility between Chlorine and Kamiokande exper- 
iments data, as well as problems for predicted berilium 
and borum neutrinos in the gallium experiments 
Recently, it was realized that by changing the solar model 
it is hardly possible to solve these problems ||. There- 
fore, it is interesting to find a solution beyond the Stan- 
dard Electroweak Model. The only known natural solu- 
tions of that kind today are the energy dependent MSW 
neutrino transitions in the Sun interior [ffol and the "just- 
so" vacuum oscillations solutions, as well as the recently 
developed hybrid solutions of MSW transitions + vacuum 
oscillations type Q . 

(b) Atmospheric neutrino anomaly: Three of the five 



underground experiments on atmospheric neutrinos have 
observed disappearance of muon neutrinos ||. This is 
in contradiction with the theoretically expected flux of 
muon neutrinos from primary cosmic rays interacting in 
the atmosphere. A successful oscillatory solution of that 
problem requires large mixing and <5m 2 of the order of 
Hr 2 eV 2 . 

(c)-Los Alamos LSND experiment claimed evidence for 
the oscillation of into v e , with a maximal probability 
of the order of 0.45 x 10~ 2 . A complementary &/„ into v e 
oscillation search, with completely different systematics 
and backgrounds, also shows a signal, which indicates the 
same favoured region of oscillation parameters . 

There exists yet another observational suggestion for 
massive neutrinos and oscillations - the dark matter prob- 
lem. Present models of structure formation in the Uni- 
verse indicate that the observed hierarchy of structures is 
reproduced best by an admixture of about 20% hot dark 
matter to the cold one fl2| |. Light neutrinos with mass 
in eV range are the only particle dark matter candidates, 
that are actually known to exist and are the most plau- 
sible candidates provided by particle physics |13|. Ac- 
tually, recent most popular hot plus cold dark matter 
models assume that two nearly degenerate massive neu- 
trinos each with mass 2.4 eV play the role of the hot dark 
matter. This small mass value is now accessible only by 
oscillations. 

However, in case we take seriously each of these exper- 
iments pointing to a neutrino anomaly and the neutrino 
oscillation solution to them, a fourth neutrino seems in- 
evitable. In the case of only three species of light neu- 
trinos with normal interactions and a see-saw hierarchy 
between the three masses, it is hardly possible to accom- 
modate all the present data simultaneously. The suc- 
cessful attempts to reconcile the LSND results with neu- 
trino oscillation solutions to the solar and atmospheric 
neutrino problems usually contain some "unnatural" fea- 
tures, like forth ultra-light sterile neutrino species, or in- 
verted neutrino mass hierarchy ]l3 ]. However, an addi- 
tional light (with mass less than 1 MeV) flavour neutrino 
is forbidden both from cosmological considerations and 
the experiments on Z decays at LEP pl|. Hence, it is 
reasonable to explore in more detail the possibility for 
an additional light sterile neutrino. Besides, GUT theo- 



ries (50(10), E 6 , etc.) |1J,|17[ and SUSY theories |8|Jl9 
predict the existence of a sterile neutrino. Moreover, 
recently models of singlet fermions, which explain the 
smallness of sterile neutrino mass and its mixing with 
the usual neutrino were proposed [fL9| . Therefore, it may 
be very useful to obtain more precise information about 
the cosmologically allowed range for the neutrino mixing 
parameters and thus present an additional independent 
test for the already discussed neutrino puzzles. More- 
over, the very small values of mass differences, which can 
be explored by the oscillations cosmological effects (like 
the ones discussed in our model) are beyond the reach of 
present and near future experiments. 

The present work is a step towards this: we suppose 
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the existence of a sterile neutrino (S , t/(2)-singlet) v s , and 
explore the cosmological effect of nonresonant neutrino 
oscillations v e «-> v s on the primordial nucleosynthesis, 
obtaining thus cosmological constraints on the neutrino 
mixing parameters. The nonresonant case in the Early 
Universe medium corresponds to the resonant case in the 
Sun, therefore, the obtained information is also of inter- 
est for the MSW solution to the solar neutrino problem. 
We discuss the special case of nonequilibrium oscillations 
between weak interacting and sterile neutrinos for small 
mass differences dm 2 , as far as the case of large Sm 2 
is already sufficiently well studied ^0|- ^7j. Oscillations 
between active and sterile neutrinos, effective before neu- 
trino freezing at 2 MeV, leading to v s thermalization be- 
fore 2 MeV have been studied there. T.e. mainly the 
equilibrium oscillations were considered with rates of os- 
cillations and neutrino weak interactions greater than the 
expansion rate. Here we discuss nonequilibrium oscilla- 
tions between electron neutrinos v e and sterile neutrinos 
v a for the case when v s do not thermalize till v e decou- 
pling at 2 MeV and oscillations become effective after 
v e decoupling. Such kind of active-sterile neutrino os- 
cillations in vacuum was first precisely studied in ]2§| ] 
using the accurate kinetic approach for the description 
of oscillating neutrinos, proposed in the pioneer work of 
Dolgov (^9|. However, the thermal background in the 
prenucleosynthesis epoch may strongly affect the propa- 
gation of neutrino pUf^l and the account of the neu- 
trino interactions with the primeval plasma is obliga- 
tory (3^]2(],|24| . The precise kinetic consideration of oscil- 
lations in a medium was provided in p4"| . It was proved 
that in case when the Universe expansion, the oscillations 
and the neutrino interactions with the medium have com- 
parable rates, their effects should be accounted for simul- 
taneously, using the exact kinetic equations for the neu- 
trino density matrix. Moreover, for the nonequilibrium 
oscillations energy distortion and asymmetry between 
neutrinos and antineutrinos may play a considerable role. 
As far as both neutrino collisions and active-sterile neu- 
trino oscillations distort the initially equilibrium active 
neutrino momentum distribution, the momentum degree 
of freedom in the description of neutrino must be ac- 
counted for. Therefore, for the case of nonequilibrium 
oscillations the evolution of neutrino ensembles should be 
studied using the exact kinetic equations for the density 
matrix of neutrinos in momentum space. This approach 
allows an exact inve stigati on of the different effects of 
neutrino oscillations ||28| , |34| , j35| : depletion of the neutrino 
number densities, the energy distortion and the genera- 
tion of asymmetry, for each separate momentum of the 
neutrino ensembles. 

In the present work we expand the original investiga- 
tion |34| for the full parameter space of the nonequilib- 
rium oscillations model for the nonresonant case. (The 
resonant case will be discussed in a following publica- 
tion.) We have provided an exact kinetic analysis of the 
neutrino evolution by a numerical integration of the ki- 
netic equations for the neutrino density matrix for each 



momentum mode. The kinetic equations are coupled non- 
linear and, therefore, an analytic solution is hardly pos- 
sible in the general case of oscillations in a medium. We 
have numerically described the evolution of the neutrino 
ensembles from the v e freezing at 2 MeV till the forma- 
tion of helium-4. 

We have calculated the production of helium-4 in a de- 
tail model of primordial nucleosynthesis, accounting for 
the direct kinetic effects of oscillations on the neutron- 
to-proton transitions. The oscillations effect on CN has 
been considered by many authors ^0|- [^Tj, |3(|, j36|- 
However, mainly the excitation of an additional de- 
gree of freedom due to oscillations (i.e. an increase of 
the effective degrees of freedom g) and the corresponding 
increase of the Universe expansion rate H ~ ^fg, leading 
to an overproduction of helium-4 was discussed. The ex- 
cluded regions for the neutrino mixing parameters were 
obtained from the requirement (based on the accordance 
between the theoretically predicted and the extracted 
from observations light elements abundances) that the 
neutrino types should be less than 3.4: N u < 3.4 |pp| - 
p4j . A successful account for the electron neutrino de- 
pletion due to oscillations was first made in |2(J and [^4j . 
In the present work we have precisely calculated the influ- 
ence of oscillations on the primordially produced helium- 
4 using the exact kinetic equations in momentum space 
for the neutron number density and the density matrix 
of neutrino, instead of their particle densities. The ac- 
curate numerical analysis of oscillations effect on helium 
production within a model of nucleosynthesis with oscil- 
lations, allowed us to account precisely for the following 
important effects of neutrino oscillations: neutrino popu- 
lation depletion, distortion of the neutrino spectrum and 
the generation of neutrino-antineutrino asymmetry. This 
enabled us to investigate the zone of very small neutrino 
mass differences up to 10 -11 eV 2 , which has not been 
reached before. As a result, we have obtained constant 
helium contours in the mass difference - mixing angle 
plane for the full range of the parameter values of our 
model. No matter what will be the preferred primordial 
helium value, favoured by future observations, it will be 
possible to obtain the excluded region of the mixing pa- 
rameters using the results of this survey. 

The paper is organized as follows. In Section II we 
present the model of nonequilibrium neutrino oscilla- 
tions. In Section III an exact analysis of the neutrino 
evolution using kinetic equations for the neutrino density 
matrix for each momentum mode is provided. The main 
effects of nonequilibrium oscillations are revealed. In Sec- 
tion IV we investigate v e into v s oscillations effect on the 
primordial production of helium using a numerical nu- 
cleosynthesis code. We discuss the influence of nonequi- 
librium neutrino oscillations, namely electron neutrino 
depletion, neutrino spectrum distortion and the genera- 
tion of neutrino-antineutrino asymmetry on the primor- 
dial yield of helium-4. The results and conclusions are 
presented in Section V. 
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II. NONEQUILIBRIUM NEUTRINO 
OSCILLATIONS - THE MODEL 

The model of nonequilibrium oscillations between weak 
interacting electron neutrinos v e and sterile neutrinos v s 
for the case when v s do not thermalize till v e decoupling 
at 2 MeV and oscillations become effective after v e decou- 
pling is described in detail in ]34] . The main assumptions 
are the following: 

• Singlet neutrinos decouple much earlier, i.e. at a 
considerably higher temperature than the active 
neutrinos do: > T^. 

This is quite a natural assumption, as far as sterile neu- 
trinos do not participate into the ordinary weak interac- 
tions. In the models predicting singlet neutrinos, the 
interactions of v s are mediated by gauge bosons with 
masses M = 0(1 TeV) [^6|jl|Q. Therefore, in later 
epochs after their decoupling, their temperature and 
number densities are considerably less than those of the 
active neutrinos due to the subsequent annihilations and 
decays of particles that have additionally heated the non- 
decoupled v e in comparison with the already decoupled 
v s . 

• We consider oscillations between v s (y s = vl) 
and the active neutrinos, according to the Majo- 
rana&Dirac (MSzD) mixing scheme [Q with mix- 
ing present just in the electron sector v$ = Uu ui, 
l = e,a:fl 

V\ = cv e + sv s 
v 2 = -sv e + cv s , 

where v s denotes the sterile electron antineutrino, c = 
cos(i9), s = sin($) and $ is the mixing angle in the elec- 
tron sector, the mass eigenstates v\ and v 2 are Majorana 
particles with masses correspondingly mi and m,%. We 
consider the nonresonant case 5m 2 = m 2 , — m\ > 0, which 
corresponds in the small mixing angle limit to a sterile 
neutrino heavier than the active one. 

In this model the clement of nonequilibrium is intro- 
duced by the presence of a small singlet neutrino density 
at 2 MeV n Vs <C n„ c , when the oscillations between v 8 
and v & become effective. In order to provide such a small 
singlet neutrino density the sterile neutrinos should have 
decoupled from the plasma sufficiently early in compari- 
son to the active ones and should have not regained their 
thermal equilibrium till 2 MeV . Therefore, as 



x The transitions between different neutrino flavours were 
proved to have negligible effect on the neutrino number den- 
sities and on primordial nucleosynthesis because of the very 
slight deviation from equilibrium in that case Tj ~ T'f (/ is 
the flavour index) §|,|g,|l). 



far as the oscillations into v e and the following noncoher- 
ent scattering off the background may lead to the ther- 
malization of v s , two more assumptions are necessary for 
the nonequilibrium case to have place: 

• Neutrino oscillations should become effective after 
the decoupling of the active neutrinos, r osc > H for 
T < 2 MeV, which is realizable for Sm 2 < 1.3 x 10" 7 
eV 2 0. 

• Sterile neutrinos should not thermalize till 2 MeV 
when oscillations become effective, i.e. the produc- 
tion rate of v s must be smaller than the expansion 
rate. 

The problem of sterile neutrino thermalization was dis- 
cussed in the pioneer work of Manohar pdj ] and in more 
recent publications [^0|- ^4|. This assumption limits 
the allowed range of oscillation parameters for our model: 
sin 2 (2i?)<5m 2 < 10" 7 eV 2 @. 

We have assumed here that electron neutrinos decou- 
ple at 2 MeV. However, the neutrino decoupling process 
is more complicated. It has been discussed in literature 
in detail |45|| . Decoupling occurs when the neutrino weak 
interaction rate r„, ~ E 2 n u (E) becomes less than the ex- 
pansion rate H ~ ^fgT 2 . Really, for electron neutrinos 
this happens at about 2 MeV. Nevertheless, due to the 
fact that weak interaction rate is greater at a higher en- 
ergy, some thermal contact between neutrinos and high 
energy plasma remains after 2 MeV, especially for the 
high energy tail of the neutrino spectrum. In case these 
high energy neutrinos begin to oscillate before their de- 
coupling, the account of this dependence of decoupling- 
time on the neutrino momentum will be essential for our 
model. Otherwise, in case these neutrinos do not start os- 
cillating before decoupling, there will be no harm consid- 
ering them decoupled earlier, as far as they preserve their 
equilibrium distribution anyway due to their extremely 
small mass. In |34j we have checked that neutrinos from 
high-energy tail start to oscillate much later than they 
decouple for the range of oscillation parameters consid- 
ered in our model. It can easily be understood from 
the fact that the oscillation rate decreases with energy 
r osc ~ Sm 2 1 E v and, therefore, neutrinos with higher en- 
ergies begin to oscillate later, namely when r osc exceeds 
the expansion rate H ~ y/gT 2 . Hence, the precise ac- 
count for the momentum dependence of the decoupling 
does not change the results of our model but unneces- 
sarily complicates the analysis and leads to an enormous 
increase of the calculation time. Therefore, in what fol- 
lows we have assumed a fixed decoupling time instead 
of considering the real decoupling period - i.e. we have 
accepted that the electron neutrinos have completely de- 
coupled at 2 MeV. 
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III. THE KINETICS OF NONEQUILIBRIUM 
NEUTRINO OSCILLATIONS 

The exact kinetic analysis of the neutrino evolution, 
discussed in this Section, though much more complicated, 
reveals some important features of nonequilibrium oscil- 
lations, that cannot be caught otherwise. As far as for the 
nonequilibrium model discussed the rates of expansion of 
the Universe, neutrino oscillations and neutrino interac- 
tions with the medium may be comparable, we have used 
kinetic equations for neutrinos accounting simultaneously 
for the participation of neutrinos into expansion, oscil- 
lations and interactions with the medium. All possible 
reactions of neutrinos with the plasma were considered, 
namely: reactions of neutrinos with the electrons, neu- 
trons and protons, neutrinos of other flavours, and the 
corresponding antiparticles, as well as self interactions of 
electron neutrinos. These equations contain all effects 
due to first order on medium-induced energy shifts, 
second order effects due to non- forward collisions, and the 
effects non-linear on the neutrino density matrices like 
neutrino refraction effects in a medium of neutrinos. In 
the case of nonequilibrium oscillations the density matrix 
of neutrinos may considerably differ from its equilibrium 
form. Then, for the correct analysis of nonequilibrium 
oscillations, it is important to work in terms of density 
matrix of neutrinos in momentum space ps|j2^ , |3^ , |35]] . 
Therefore, we have provided a proper kinetic analysis 
of the neutrino evolution using kinetic equations for the 
neutrino density matrix for each momentum mode. 

Hence, the kinetic equations for the density matrix of 
the nonequilibrium oscillating neutrinos in the primeval 
plasma of the Universe in the epoch previous to nucle- 
osynthesis, i.e. consisting of photons, neutrinos, elec- 
trons, nucleons, and the corresponding antiparticles, 
have the form: 



dpit) 
dt 



Hp 



d P (t) 
dp 



■i[H ,p(t)} + i[H mt ,p(t)} 



{Hj nt ) 
(1) 



where p is the momentum of electron neutrino and p is 
the density matrix of the massive Majorana neutrinos in 
momentum space. 

The first term in the equation describes the effect of 
expansion, the second is responsible for oscillations, the 



2 When neutrinos are in equilibrium their density matrix has 
its equilibrium form, namely pij = 8% 3 exp(/x/T — E/T), so 
that one can work with particle densities instead of p. In 
an equilibrium background, the introduction of oscillations 
slightly shifts p from its diagonal form, due to the extreme 
smallness of the neutrino mass in comparison with the char- 
acteristic temperatures and to the fact that equilibrium dis- 
tribution of massless particles is not changed by the expan- 
sion Ei. 



third accounts for forward neutrino scattering off the 
medium and the last one accounts for second order inter- 
action effects of neutrinos with the medium. TL Q is the 
free neutrino Hamiltonian: 



H =( & 







^p 1 + ml 



while TLi n t = a V is the interaction Hamiltonian, where 
aiij = U* e U je , V = Gf (+L - Q/Myy), and in the inter- 
action basis plays the role of an induced squared mass 
for electron neutrinos: 



n-f LR — 

ri int — 



V 




Hence, V is the time varying (due to the Universe cool- 
ing) effective potential, induced by the interactions of 
neutrino with the medium through which it propagates. 
Since v s does not interact with the medium it has no 
self-energy correction, i.e. V s = 0. 

The first 'local' term in V accounts for charged- and 
neutral-current tree-level interactions of v e with medium 
protons, neutrons, electrons and positrons, neutrinos and 
antincutrinos. It is proportional to the fermion asymme- 
try of the plasma L — Lj, which is usually taken to 
be of the order of the baryon one i.e. 10~ 10 (i.e. B — L 
conservation is assumed). 



L 



f 



N f - N f 



rpd 



N B -N £ 
N~ 



r^6 



f3T 3 . 



The second 'nonlocal' term in V arises as an W/Z prop- 
agator effect, Q ~ E v T 4 j33| , pfH . For the early Universe 
conditions both terms must be accounted for because al- 
though the second term is of the second power of Gp , the 
first term is proportional besides to the first power of Gf, 
also to the small value of the fermion asymmetry. More- 
over, the two terms have different temperature depen- 
dence and an interesting interplay between them during 
the cooling of the Universe is observed. At high temper- 
ature the nonlocal term dominates, while with cooling of 
the Universe in the process of expansion the local one 
becomes more important. 

The last term in the Eq. (||) describes the weak interac- 
tions of neutrinos with the medium. For example, for the 
weak reactions of neutrinos with electrons and positrons 
e + e~ <-» Vii/j, e^v.\ — ► e'^v' it has the form 



/ dQ(£,e+,e-) 
+ Jdft(e-y,e'-) 



n e ~n e +AA^ - \ {p, A ] pA} + 
> e+ Cp'tf -i{CtC, p},n e+ 



where n stands for the number density of the interacting 
particles, 



<m(i,j,k) 



(2^) 4 r d 3 pi 

2E V J (2tt) ;1 2Ei (2tt) 3 1Ej (27r) 3 " 2E k 



Pi 



y-8 A (p v +pi -pj -pk) 
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is a phase space factor, A is the amplitude of the pro- 
cess e + e~ — > Viiij, B is the amplitude of the process 
eTvj — > e'~i/ and C is the amplitude of the process 
They are expressed through the known 



amplitudes A e (e + e 
C e (e + v e — * e + v e ): 

v4 — j^p. 



f e i? e ), ^e(e f e — * e i/ e ) and 



B = a B e , C — a C e 



An analogues equations hold for the antineutrino den- 
sity matrix, the only difference being in the sign of the 
lepton asymmetry: Lf is replaced by — Lf. Medium 
terms depend on neutrino density, thus introducing a 
nonlinear feedback mechanism. Neutrino and antineu- 
trino ensembles evolve differently as far as the back- 
ground is not CP symmetric. Oscillations may change 
neutrino-antineutrino asymmetry and it in turn affects 
oscillations. The evolution of neutrino and antineutrino 
ensembles is coupled and hence, it must be considered 
simultaneously. 

We have analysed the evolution of the neutrino density 
matrix for the case when oscillations become noticeable 
after electron neutrinos decoupling, i.e. after 2 MeV. 
Then the last term in the kinetic equation can be ne- 
glected. So, the equation ([!]) results into a set of coupled 
nonlinear integro-differential equations with time depen- 
dent coefficients for the components of the density matrix 
of neutrino, ft is convenient instead of d/dt to use d/d/j,, 
where /i 2 = ^/f67r 3 g/45 (S 2 /Mpi) t, and S = m n — m p . 




Then from eq. (Q) we obtain: 

(Ax 
p'22 
p'12 

\P'21 



where prime denotes d/d/j, and M = 5m 2 /(2E u ) + (s 2 - 
c 2 )V. 

Analytical solution is not possible without drastic as- 
sumptions and, therefore, we have numerically explored 
the problem ^| using the Simpson method for integration 
and the fourth order Runge-Kutta algorithm for the so- 
lution of the differential equations. 

The neutrino kinetics down to 2 MeV does not differ 
from the standard case, i.e. electron neutrinos maintain 
their equilibrium distribution, while sterile neutrinos are 
absent. So, the initial condition for the neutrino ensem- 
bles in the interaction basis can be assumed of the form: 



where <<z = exp(-E u /T)/(l + exp(-E v /T)). 

We have analyzed the evolution of nonequilibrium os- 
cillating neutrinos by numerically integrating the kinetic 
equations (§) for the period after the electron neutrino 
decoupling till the freeze out of the neutron-proton ra- 
tio (n/p-ratio), i.e. for the temperature interval [0.3,2.0] 
MeV. The oscillation parameters range studied is 8m 2 £ 
[10- n , 10- J ] eV 2 and i9 £ [0, tt/4]. 

The distributions of electrons and positrons were taken 
the equilibrium ones. Really, due to the enormous rates 
of the electromagnetic reactions of these particles the de- 
viations from equilibrium are negligible. We have also 
neglected the distortion of the neutrino spectra due to 
residual interactions between the electromagnetic and 
neutrino components of the plasma after 2 MeV. This 
distortion was accurately studied in JfT| , where it was 
shown that the relative corrections to v e density is less 
than 1 % and the effect on the primordial helium abun- 
dance is negligible. 

The neutron and proton number densities, used in the 
kinetic equations for neutrinos, were substituted from the 
numnerical calculations in CN code accounting for neu- 
trino oscillations. I.e. we have simultaneously solved 
the equations governing the evolution of neutrino en- 
sembles and those describing the evolution of the nu- 
cleons (see the next section). The baryon asymmetry /?, 
parametrized as the ratio of the baryon number density 
to the photon number density, was taken to be 3 x 10 -10 . 

Three main effects of neutrino nonequilibrium oscil- 
lations were revealed and precisely studied, namely elec- 
tron neutrino depletion, neutrino energy spectrum distor- 
tion and the generation of asymmetry between neutrinos 
and their antiparticles: 

(a) Depletion of v e population due to oscillations: As 
far as oscillations become effective when the number den- 
sities of v e are much greater than those of v s , N Ue ^> N Vs , 
the oscillations tend to reestablish the statistical equilib- 
rium between different oscillating species. As a result 
N Va decreases in comparison to its standard equilibrium 
value due to oscillations in favour of sterile neutrinos. 
The effect of depletion may be very strong (up to 50%) 
for relatively great Sm 2 and maximal mixing. This re- 
sult of our study is in accordance with other publications 
concerning depletion of electron neutrino population due 
to oscillations, like p^j, however, we have provided more 
precise account for this effect due to the accurate kinetic 
approach used. 

In Figs. 1 the evolution of neutrino number densities 
is plotted. In Fig. la the curves represent the evolu- 
tion of the electron neutrino number density in the dis- 
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3 For the case of vacuum neutrino oscillations this equation 
was analytically solved and the evolution of density matrix 
was given explicitly in Ref. Eq]. 



4 Note, that while neutrinos are in thermal equilibrium with 
the plasma no dilution of their number density is expected as 
far as it is kept the equilibrium one due to the annihilations 
of the medium electrons and positrons. 
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cussed model with a fixed mass difference 5m 2 = 10~ 8 
eV 2 and for different mixings. The numerical analysis 
showed that for small mixing, sin 2 (2i9) < 0.01, the re- 
sults do not differ from the standard case, i.e. then oscil- 
lations may be neglected. In Fig. lb the evolution of the 
electron neutrino number density is shown for a nearly 
maximum mixing, sin 2 (2$) = 0.98, and different squared 
mass differences. Our analysis has proved, that for mass 
differences <5m 2 < 10~ n eV 2 , the effect of oscillations is 
negligible for any <&. 

In case of oscillations effective after the neutrino freeze 
out, electron neutrinos are not in thermal contact with 
the plasma and, therefore, the electron neutrino state, de- 
pleted due to oscillations into steriles, cannot be refilled 
by electron-positron annihilations. That irreversible de- 
pletion of v e population exactly equals the increase of v a 
one (see Fig. lc). The number of the effective degrees of 
freedom do not change due to oscillations in that case, 
as far as the electron neutrino together with the corre- 
sponding sterile one contribute to the energy density of 
the Universe as one neutrino unit, even in case when the 
steriles are brought into chemical equilibrium with i/ e . 
This fact was first noted in ]23]].p| 

(b) Distortion of the energy distribution of neutrinos: 
The effect was first discussed in Q for the case of flavour 
neutrino oscillations. However, as far as the energy dis- 
tortion for that case was shown to be negligible p9| , p| , 
it was not paid the necessary attention it deserved. The 
distortion of the neutrino spectrum was not discussed 
in publications concerning active-sterile neutrino oscilla- 
tions, and was thought to be negligible. In [^8| it was first 
shown that for the case of v e <-> v s vacuum oscillations 
this effect is considerable and may even exceed that of 
an additional neutrino species. In |34j we have discussed 
this effect for the general case of neutrino oscillations in a 
medium. The evolution of the distortion is the following: 
Different momentum neutrinos begin to oscillate at dif- 
ferent temperatures and with different amplitudes. First 
the low energy part of the spectrum is distorted, and 
later on this distortion concerns neutrinos with higher 
and higher energies. This behaviour is natural, as far 
as neutrino oscillations affect first low energy neutrinos, 
r osc ~ dm 2 /E v . The Figs. 2a, 2b, 2c and 2d snapshot 
the evolution of the energy spectrum distortion of active 
neutrinos x 2 pll(x), where x = E„/T, for maximal mix- 
ing and Sm 2 = 10~ 85 eV 2 , at different temperatures: 
T = 1 MeV (a), T = 0.7 MeV (b), T = 0.5 MeV (c), 
T = 0.3 MeV (d). As can be seen from the figures, the 
distortion down to temperatures of 1 MeV is not signif- 
icant as far as oscillations are not very effective and/or 



5 Note the essential difference from the case of electron neu- 
trinos in thermal equilibrium, when the oscillations into sterile 
neutrinos bring an additional degree of freedom into thermal 
contact. 



the weak residual interactions with the background still 
can compensate for the difference. However, for lower 
temperatures the distortion increases and at 0.5 MeV is 
strongly expressed. Its proper account is important for 
the correct determination of oscillations role in the kinet- 
ics of n-p transitions during the freeze out of nucleons at 
about 0.3 MeV. 

Our analysis has shown that the account for the 
nonequilibrium distribution by shifting the effective tem- 
perature and assuming the neutrino spectrum of equilib- 
rium form, often used in literature (see for example pG| ) , 
may give misleading results for the case 8m 2 < 10 -7 eV 2 . 
The effect cannot be absorbed merely in shifting the effec- 
tive temperature and assuming equilibrium distributions. 
For larger neutrino mass differences oscillations are fast 
enough and the naive account is more acceptable, pro- 
vided that v e have not decoupled. 

(c) The generation of asymmetry between v e and 
their antiparticles: The problem of asymmetry genera- 
tion in different contexts was discussed by several au- 
thors. The possibility of an asymmetry generation due 
to CP-violating flavour oscillations was first proposed in 
Ref. ]3§| ] . Later estimations of an asymmetry due to CP- 
violating MSW resonant oscillations were provided [ p9| . 
The problem of asymmetry was considered in connection 
with the exploration of the neutrino propagation in the 
early Universe CP-odd plasma also in pQ] - H] and this 
type of asymmetry was shown to be negligible. Recently 
it was realized in |40 27]], that asymmetry can grow to a 
considerable values for the case of great mass differences, 
Sm 2 > 10~ 5 eV 2 . The effect of asymmetry for small 
mass differences Sm 2 < 10~ 7 eV 2 on primordial produc- 
tion of helium was also proved to be important for the 
case of resonant neutrino oscillations (3^] . Our approach 
allows precise description of the asymmetry evolution, as 
far as working with the self consistent kinetic equations 
for neutrinos in momentum space enables us to calcu- 
late the behaviour of the asymmetry at each momentum. 
This is important particularly when the distortion of the 
neutrino spectrum is considerable. 

In the present work we have explored accurately the 
effect of the asymmetry in the nonresonant case for all 
mixing angles and for small mass differences 5m < 10~ 7 
eV 2 . Our analysis showed that when the lepton asym- 
metry is accepted initially equal to the baryon one, (as is 
usually assumed for the popular L — B conserving mod- 
els), the effect of the asymmetry is small for all the dis- 
cussed parameters range. And although the asymmetry 
is not wiped out by the coupled oscillations, as stated by 
some authors nonresonant neutrino oscillations 

really cannot generate large neutrino-antineutrino asym- 
metry in the early Universe. This result is in accordance 
with the conclusions concerning asymmetry evolution in 
34 1 . We have also checked that the neutrino asym- 



metry even in the case of initial neutrino asymmetry by 
two orders of magnitude higher does not have significant 
effect on the cosmologically produced 4 He. Therefore, 
for such small initial values of the lepton asymmetry, the 
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neutrino asymmetry should be better neglected when cal- 
culating primordial element production for the sake of 
computational time. Mind, however, that for higher val- 
ues of the initial asymmetry the effect could be signifi- 
cant, and should be studied in detail. The asymmetry 
evolution and its effect on He-4 production for unusual 
high initial values of the lepton asymmetry will be stud- 
ied elsewhere f46| . 

In conclusion, our numerical analysis showed that the 
nonequilibrium oscillations can considerably deplete the 
number densities of electron neutrinos (antincutrinos) 
and distort their energy spectrum. 



IV. NUCLEOSYNTHESIS WITH 
NONEQUILIBRIUM OSCILLATING NEUTRINOS 

As an illustration of the importance of these effects, 
and hence of the proposed approach to the analysis of 
nonequilibrium neutrino oscillations, we discuss their in- 
fluence on the primordial production of 4 He. The effect 
of oscillations on nucleosynthesis has been discussed in 



39|. A detail ki- 
ium for the case 



numerous publications []20||- J28| , |3l 
netic calculation of primordial yield of he 
of the nonequilibrium oscillations in vacuum was made 
in p^j and the proper consideration accounting for the 
neutrino forward scattering processes off the background 
particles was done in |34[ for some neutrino mixing pa- 
rameters. [] In the present work we calculate precisely 
the influence of oscillations on the production of He-4 
within a detail numerical CN model with nonresonant 
nonequilibrium neutrino oscillations. The analysis of ]34| ] 
is expanded for the full space of the mixing parameters 
values. 

Working with exact kinetic equations for the nucleon 
number densities and neutrino density matrix in momen- 
tum space, enables us to analyze the direct influence of 
oscillations onto the kinetics of the neutron-to-proton 
transfers and to account precisely for the neutrino de- 
pletion, neutrino energy distortion and the generation of 
asymmetry due to oscillations. 

Primordial element abundances depend primarily on 
the ncutron-to-proton ratio at the weak freeze out 
((n/p) /-ratio) of the reactions interconverting neutrons 
and protons : n + v e <-* p + e and n + e + <-> p + i7 e . 
The freeze out occurs when due to the decrease of tem- 
perature with Universe expansion these weak interaction 



6 Calculations of helium production within the full Big 
Bang Nucleosynthesis code with oscillations were provided 
also in (26) , however, there the momentum degree of freedom 
of neutrino was not considered and a simplifying account of 
the nonequilibrium was used - by merely shifting the neutrino 
effective temperature and working in terms of equilibrium par- 
ticle densities. 



rates r,„ ~ E 2 n v become comparable and less than the 
expansion rate H ~ yfg T 2 . Hence, the (n/p)/-ratio de- 
pends on the effective relativistic degrees of freedom g 
(through the expansion rate) and the neutrino number 
densities and neutrino energy distribution (through the 
weak rates). Therefore, we calculate accurately the evo- 
lution of neutron number density till its freeze-out. Fur- 
ther evolution is due to the neutron decays n — > p+e + i7 e 
that proceed till the effective synthesis of deuterium be- 
gins. As far as the expansion rate exceeds considerably 
the decay rate for the characteristic period before the 
freeze out, decays are not essential. Therefore, we have 
accounted for them adiabatically. 

The master equation, describing the evolution of the 
neutron number density in momentum space n n for the 
case of oscillating neutrinos v e «-> v s , reads: 



(dUn/dt) = Hp n (drin/dp n ) + 

+ / dil(e-,p,v)\A(e-p-> vn)\ 2 

x \n K -n p (l - pll) - n n p LL (l 



Jdfl(e+,p,D)\A(e^ 
x [n e +n„(l - pll) 



pv)\ 2 
- n p p L L(l ~ n e +)] 



■)] (3) 



The first term on the right-hand side describes the ef- 
fect of expansion while the next ones - the processes 
e~ + p «-> n + v e and p + v e <-> e + + n, directly influ- 
encing the nucleon density. It differs from the standard 
scenario one only by the substitution of pll and pll 
instead of nf, q = [1 — exp(7S„ /7 1 )] -1 . The neutrino and 
antineutrino density matrices differ Pll 7^ Pll, contrary 
to the standard model, as a result of the different reac- 
tions with the CP-odd plasma of the prenucleosynthesis 
epoch. We have accounted for the final state Pauli block- 
ing for neutrinos and electrons. 

Particle number densities per unit volume are ex- 
pressed as N — (2tt)~ 3 J d 3 p n(p). Performing the in- 
tegration on the right-hand side of the equation also one 
gets the final equations for the time evolution of the neu- 
tron number density: 

(dNjdt) = -3HN n + G% g v T 5 x 

7T 

x [N p J [I- pll(x)} 1 + e _ x _ y fix, y)dx 

,00 ! 
- N n pll{x) 1 - e _ x _ y f{x,y)dx 



+ N P pll{x) f(x,-y)dx 

J(i+Ov 1 + e x+y 



(i+Ow 



[1 - pll{x)} 



p -x+y 



1 + e~ x +y 



f(x,-y)dx] (4) 



where f(x, y) = x 2 (x + y) yj(x + y) 2 + C, 2 y 2 and y = (5 + 
m e )/T, C = m e /S, 5 = m n - m p . 

The first term on the right-hand side describes the di- 
lution effect of expansion, the next describe the weak 
processes, as pointed above. We have numerically inte- 
grated this equation for the temperature range of interest 
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T £ [0.3, 2.0] MeV for the full range of oscillation parame- 
ters of our model. The value of pll{x) at each integration 
step was taken from the simultaneously performed inte- 
gration of the set of equations (||), i.e. the evolution of 
neutrino and the nuclcons was followed self consistently. 
The initial values at T = 2 MeV for the neutron, proton 
and electron number densities are their equilibrium val- 
ues. Although the electron mass is comparable with the 
temperature in the discussed temperature range, the de- 
viation of the electron density from its equilibrium value 
is negligible due to the enormous rate of the reactions 
with the plasma photons [^9| . The parameters values of 
the CN model, adopted in our calculations, are the fol- 
lowing: the mean neutron lifetime is t = 887 sec, which 
corresponds to the present weighted average value jl7), 
the effective number of relativistic flavour types of neu- 
trinos during the nucleosynthesis epoch N v is assumed 
equal to the standard value 3. This is a natural choice 
as far as it is in good agreement both with the CN argu- 
ments |Q [] and with the precision measurements of the 
Z decay width at LEP JOf. 

V. RESULTS AND CONCLUSIONS 

The results of the numerical integration are illustrated 
in Fig. 3. As it can be seen from the figure the kinetic 
effects (neutrino population depletion and distortion of 
neutrino spectrum) due to oscillations play an important 
role and lead to a considerable overproduction of helium. 

Qualitatively the effect of oscillations on helium pro- 
duction can be described as follows: 

The depletion of the electron neutrino number densi- 
ties due to oscillations into sterile ones strongly affects 
the n <-> p reactions rates. It leads to an effective de- 
crease in the processes rates, and hence to an increase 
of the freezing temperature of the n/p-ratio and the cor- 
responding overproduction of the primordially produced 

The effect of the distortion of the energy distribution 
of neutrinos has two aspects. On one hand an average 
decrease of the energy of active neutrinos leads to a de- 
crease of the weak reactions rate, T w ~ E 2 and subse- 
quently to an increase in the freezing temperature and 
the produced helium. On the other hand, there exists 
an energy threshold for the reaction v e + p — > n + e + . 
And in case when, due to oscillations, the energy of the 
relatively greater part of neutrinos becomes smaller than 
that threshold the n/p- freezing ratio decreases leading 
to a corresponding decrease of the primordially produced 
helium-4 [M. The numerical analysis showed that the 



7 However mind also the possibilities for somewhat relaxation 
of that kind of bound in modifications of the CN model with 
decaying particles as in |4S|,pO| . 



latter effect is less noticeable compared with the former 
ones. 

The asymmetry calculations showed a slight predom- 
inance of neutrinos over antineutrinos, not leading to 
a noticeable effect on the production of helium in case 
the lepton asymmetry is accepted initially equal to the 
baryon one. So, the effect of asymmetry is proved to 
be negligible for all the discussed parameter range, i.e. 
for any d and for 5m 2 < 10 -7 eV 2 . We have partially 
(not for the full range of model parameters) investigated 
the problem for higher than the baryon one initial lep- 
ton asymmetry. The preliminary results point that even 
lepton asymmetry initially by two orders of magnitude 
higher does not have noticeable effect on the cosmologi- 
cally produced 4 He. Higher than those lepton asymme- 
tries, however, should be accounted for properly even in 
the nonresonant case. 

Thus, the total result of nonequilibrium neutrino oscil- 
lations is an overproduction of helium in comparison to 
the standard value. 

In Fig. 4 the dependence of the frozen neutron number 
density relative to nucleons X n = N n / (N p + N n ) on the 
mixing angle for different fixed Sm 2 is illustrated. The 
dependence of the frozen neutron number density rela- 
tive to nucleons X n — N n /(N P + N n ) on the Sm 2 for 
fixed different mixing angles, is presented in Fig. 5. The 
effect of oscillations is maximal at maximal mixing for 
the nonresonant case of neutrino oscillations. As it can 
be seen from the figures, it becomes almost negligible 
(less than 1%) for mixings as small as 0.1 for any Sm 2 of 
the discussed range of our model. The value of the frozen 
n/p-ratio is a smoothly increasing function of the mass 
difference. Our analysis shows that the effect of oscilla- 
tion for 8m 2 smaller than 10 -10 eV 2 even for maximal 
mixing is smaller than 1%. The nonresonant oscillations 
with 8m 2 < 10 -11 eV 2 do not have any observable effect 
on the primordial production of elements, i.e. the re- 
sults coincide with the standard model values with great 
accuracy. 

From the numerical integration for different oscilla- 
tion parameters we have obtained the primordial helium 
yield Y p (5m 2 , which is illustrated by the surface in 
Fig. 6. Some of the constant helium contours calculated 
in the discussed model of cosmological nucleosynthesis 
with nonresonant neutrino oscillations on the Sm 2 — d 
plane are presented in Fig. 7. 

On the basis of these results, requiring an agreement 
between the theoretically predicted and the observa- 
tional values of helium, it is possible to obtain cosmo- 
logical constraints on the neutrino mixing parameters. 
At present the primordial helium values extracted from 
observations differ considerably: for example some au- 
thors believe that the systematic errors have already 
been reduced to about the same level as the statistical 
one and obtain the bounds for the primordial helium: 
Y p ( A He) = 0.232 ±0.003 @, while others argue that un- 
derestimation of the systematic errors, such as errors in 
helium emissivities, inadequatisies in the radiative trans- 







fer model used, corrections for underlying stellar absorp- 
tion and fluorescent enhancement in the He I lines, cor- 
rections for neutral helium, may be significant and their 
account may raise the upper bound on Y p as high as 
0.26 J52j . Thus besides the widely adopted "classical" 
bound Y p < 0.24 |33| it is reasonable to have in mind 
the more "reliable" upper bound to the primordial he- 
lium abundance Y v < 0.25 and even the extreme value 
as high as 0.26 154]. Therefore, we considered it useful 
to provide the precise calculations for helium contours 
up to 0.26. So, whatever the primordial abundance of 
4 i?e will be found to be in future (within this extreme 
range) the results of our calculations may provide the 
corresponding bound on mixing parameters of neutrino 
for the case of nonresonant active-sterile oscillations with 
small mass differences. Assuming the conventional obser- 
vational bound on primordial 4 H e 0.24 the cosmologically 
excluded region for the oscillation parameters is shown on 
the plane sin 2 (2-i9) - 5m 2 in Fig. 7. It is situated to the 
right of the Y p — 0.245 curve, which gives 5% overpro- 
duction of helium in comparison with the accepted 0.24 
observational value. 

The curves, corresponding to helium abundance Y p — 
0.24, obtained in the present work, and in previous works, 
analyzing the nonresonant active-sterile neutrino oscilla- 
tions, are plotted in Fig. 8. In |^(| and the authors 
estimated the effect of excitement of an additional de- 
gree of freedom due to oscillations, and the correspond- 
ing increase of the Universe expansion rate, leading to 
an overproduction of helium-4. The excluded regions for 
the neutrino mixing parameters were obtained from the 
requirement that the neutrino types should be less than 
3.4: N v < 3.4. In these works the depletion effect was 
considered. The asymmetry was neglected and the dis- 
tortion of the neutrino spectrum was not studied as far 
as the kinetic equations for neutrino mean number densi- 
ties were considered. Our results are in good accordance 
with the estimations in [ pof and the numerical analysis 
in p3|] , who have made very successful account for one of 
the discussed effects of nonequilibrium oscillations - the 
neutrino population depletion. The results of j2(|, as can 
be seen from the Fig. 8, differ more both from the ones of 
the previously cited works and from our results. Proba- 
bly the account for nonequilibrium oscillations merely by 
shifting the effective neutrino temperature, as assumed 
there is not acceptable for a large range of model param- 
eters. 

As can be seen from the curves, for large mixing an- 
gles, we exclude Sm 2 > 10~ 9 eV 2 , which is almost an 
order of magnitude stronger constraint than the previ- 
ously existing. This more stringent constraints obtained 
in our work for the region of great mixing angles and 
small mass differences is due to the more accurate ki- 
netic approach we have used and to the precise account 
of neutrino depletion, energy distortion and asymmetry 
due to oscillations. 

As far as we already have at our disposal some impres- 
sive indications for neutrino oscillations, it is interesting 



to compare our results also with the range of parame- 
ters which could eventually explain the observed neutrino 
anomalies: 

The vacuum oscillation interpretation of the solar neu- 
trino problem requires extremely small mass differences 
squared, less and of the order of 10~ 10 eV 2 . It is safely 
lower than the excluded region, obtained in our work, 
and is, therefore, allowed from CN considerations. The 
MWS small mixing angle nonadiabatic solution (see for 
example Krastev, Liu and Petcov in jll]]) is out of the 
reach of our model. However, as we are in a good ac- 
cordance with the results of active-sterile neutrino oscil- 
lation models with higher mass differences, it is obvious 
that a natural extrapolation of our excluded zone towards 
higher mass differences will rule out partially the possible 
solution range for large mixing angles. 

Our pattern of neutrino mixing is compatible with 
models of degenerate neutrino masses of the order of 2.4 
eV, necessary for the successful modelling of the struc- 
ture formation of the Universe in Hot plus Cold Dark 
Matter Models 0. 

As a conclusion, we would like to outline the main 
achievements of this work: In a model of nonequilibrium 
nonresonant active-sterile oscillation, we had studied the 
effect of oscillations on the evolution of the neutrino num- 
ber densities, neutrino spectrum distortion and neutrino- 
antincutrino asymmetry. We have used kinetic equations 
for the density matrix of neutrinos in momentum space, 
accounting simultaneously for expansion, oscillations and 
interactions with the medium. This approach enabled 
us to describe precisely the behaviour of neutrino en- 
sembles in the Early Universe in the period of interest 
for CN. The analysis was provided for small mass differ- 
ences. We have shown that the energy distortion may 
be significant, while the asymmetry in case it is initially 
(i.e. before oscillations become effective) of the order of 
the baryon one, may be neglected. 

Next, we have made a precise survey of the influence 
of the discussed type of oscillations on the cosmological 
production of helium-4. We have calculated the evolu- 
tion of the corresponding neutron-to-proton ratio from 
the time of freeze out of neutrinos at 2 MeV till the 
effective freeze out of nucleons at 0.3 MeV for the full 
range of model parameters. As a result we have obtained 
the dependence Y p (5m 2 , i?) and constant helium contours 
on the Sm 2 — $ plane. Requiring an agreement between 
the observational and the theoretically predicted primor- 
dial helium abundances, we have calculated accurately 
the excluded regions for the neutrino mixing parameters, 
for different assumptions about the preferred primordial 
value of helium. 
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Figure Captions 



Figure la: The curves represent the calculated evolution of the electron neutrino number density in the discussed 
model of active-sterile neutrino oscillations with a mass difference Sm 2 = 10~ 8 eV 2 and different mixing, parametrized 
by sin 2 (2i9), namely: 1, 10" 01 , 1CT - 1 and 0.1. 

Figure lb: The curves show the evolution of the electron neutrino number density in the discussed model of 
nonresonant active-sterile neutrino oscillations for a nearly maximum mixing, sin 2 (2$) = 0.98, and different squared 
mass differences Sm 2 , namely 10~ 7 , 10~ 8 , 10~ 9 and 10~ 10 in cV 2 . 

Figure lc: The curves show the evolution of the electron neutrino number density (the solid curve) and the sterile 
neutrino number density (the dashed curve) in the case of the nonresonant active-sterile neutrino oscillations for a 
maximal mixing and Sm 2 = 10~ 8 cV 2 . The reduction of the active neutrino population is exactly counterbalanced by 
a corresponding increase in the sterile neutrino population. 

Figure 2: The figures illustrate the evolution of the energy spectrum distortion of active neutrinos x 2 pll{x), where 
x = E„/T, for the case of nonresonant v e -v s oscillations with a maximal mixing and Sm 2 = 10 -8 ' 5 eV 2 , at different 
temperatures: T = 1 MeV (a), T = 0.7 MeV (b), T = 0.5 MeV (c), T = 0.3 MeV (d). 

Figure 3: The evolution of the neutron number density relative to nucleons X n (t) = N n (t)/(N p + N n ) for the case 
of nonresonant oscillations with maximal mixing and different 5m 2 is shown. For comparison the standard model 
curve is plotted also. 



Figure 4: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X n = 
N n /(N P + N n ) on the mixing angle for different Sm 2 . 

Figure 5: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X n = 
N n /(N p + N n ) on the mass difference for different mixing angles. 



Figure 6: The dependence of the primordially produced helium on the oscillation parameters is represented by the 
surface Y p (Sm 2 , $) . 

Figure 7: On the Sm 2 — i9 plane some of the constant helium contours calculated in the discussed model of 
cosmological nucleosynthesis with nonresonant neutrino oscillations are shown. 

Figure 8: The curves, corresponding to helium abundance Y p — 0.24, obtained in the present work and in previous 
works, analyzing the nonresonant active-sterile neutrino oscillations, are plotted on the Sm 2 — •& plane. 




Figure la: The curves represent the calculated evolution of the electron neutrino number density in the discussed 
model of active-sterile neutrino oscillations with a mass difference Sm 2 = 10~ 8 cV 2 and different mixing, parametrized 
by sin 2 (2tf), namely: 1, 10" 001 , 1CT 01 and 0.1. 




Figure lb: The curves show the evolution of the electron neutrino number density in the discussed model of 
nonresonant active-sterile neutrino oscillations for a nearly maximum mixing, sin 2 (2$) = 0.98, and different squared 
mass differences 8m 2 , namely 10~ 7 , 1CP 8 , 10~ 9 and 1CP 10 in cV 2 . 
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Figure lc: The curves show the evolution of the electron neutrino number density (the solid curve) and the sterile 
neutrino number density (the dashed curve) in the case of the nonresonant active-sterile neutrino oscillations for a 
maximal mixing and 5m 2 = 1CP 8 cV 2 . The reduction of the active neutrino population is exactly counterbalanced by 
a corresponding increase in the sterile neutrino population. 
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Figure 2: The figures illustrate the evolution of the energy spectrum distortion of active neutrinos x 2 Pll(x), where 
x = E v /T, for the case of nonresonant v e -v s oscillations with a maximal mixing and 6m 2 = 10~ 8 5 eV 2 , at different 
temperatures: T = 1 MeV (a), T = 0.7 MeV (b), T = 0.5 MeV (c), T = 0.3 MeV (d). 
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Figure 3: The evolution of the neutron number density relative to nucleons X n {t) — N n {t)/{N p + N n ) for the case 
of nonresonant oscillations with maximal mixing and different <5m 2 is shown. For comparison the standard model 
curve is plotted also. 




Figure 4: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X n = 
N n /(N p + N n ) on the mixing angle for different 5m 2 . 




Figure 5: The figure illustrates the dependence of the frozen neutron number density relative to nucleons X, 
N n /(N p + N n ) on the mass difference for different mixing angles. 




Figure 7: On the 5m 2 — i9 plane some of the constant helium contours calculated in the discussed model of 
cosmological nucleosynthesis with nonresonant neutrino oscillations are shown. 




Figure 8: The curves, corresponding to helium abundance Y p = 0.24, obtained in the present work and in previous 
works, analyzing the nonresonant active-sterile neutrino oscillations, are plotted on the 5m 2 — •& plane. 



